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VORTEX FLOWMETER 

CROSS-REFERENCE TO RELATED APPLICATIONS 

This application claims priority from U.S. Provisional Application No. 60/261,196, 
filed January 16, 2001, and titled VORTEX FILTER, which is incorporated by reference. 

5 TECHNICAL FIELD 

This invention generally relates to digital and analog flowmeters. 

BACKGROUND 

Flowmeters provide information about materials being transferred through a conduit. 
For example, a flowmeter may provide a direct indication of the velocity and volume of the 
10 material being transferred through the conduit. Mass flowmeters measure the mass of the 
material, and density flowmeters measure the density of the material. 

One type of flowmeter is a vortex-shedding meter, or vortex flowmeter. Vortex 
flowmeters operate using techniques based on a vortex shedding phenomenon known as the 
I Karman effect. In the vortex flowmeter, a bluff body is placed in the conduit (e.g., a pipe) in 

I 15 the path of the material passing through the conduit. The bluff body sheds Karman vortices 
: in its wake. These vortices cause fluctuating pressure variations that then are amplified and 

* processed. The frequency of the Karman vortices is directly proportional to the flow rate of 

the material, as are the pressure variations occurring during vortex generation. 

To produce a flow measurement, the flowmeter first determines the vortex shedding 
20 frequency, oh in rad/sec. The flowmeter also converts the reading into a flow,/ in units of 
1/sec, through use off=K m Q)v, where K m is the Appropriate meter factor. This meter factor K m 
(or K-factor) relates the frequency of generated vortices to the flow rate. 



Conventional vortex flowmeters can accurately measure the vortex shedding 
frequency down to a threshold low flow rate of material passing through the conduit. Below 
25 this threshold, accuracy of conventional measurement becomes unreliable. This primarily 
results because the amplitude of the vortex shedding pressure is much smaller at low flow 
rates, such that the signal is much noisier due to harmonic components caused by pressure 
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pulsations from, for example, the impeller blades of the pump(s). Therefore, at low flow 
rates, the decreased amplitude of the vortex shedding pressure and the increased effect of 
noise on the measurement system lead to greater uncertainties in the flow rate measurement 
and/or slower determinations of flow rate. 



5 SUMMARY 

Techniques are provided for improving the accuracy of flow rate measurement over a 
greater range of flow rates. In particular, the techniques permit flow meters to measure low 
flow rates as well as medium and high flow rates with greater accuracy and speed. 

In one implementation, a process variable transmitter includes a vortex flowmeter 
10 using multiple phase-locked loops (PLL). In particular, the flowmeter may include a first 
PLL having a first bandwidth producing a first output signal, and a second PLL having a 
second bandwidth narrower than the first bandwidth of the first PLL. The first and second 
PLLs are operable to lock into a frequency of an input signal and produce first and second 
output signals, respectively. The second PLL is operable to lock into the frequency of the 
1 5 input signal with greater accuracy and greater immunity to noise than the first PLL. 

Implementations may include one or more of the following features. For example, 
the process variable transmitter may include a switch operable to switch an output signal of 
the process variable transmitter between the first output signal and the second output signal. 
The second PLL may generate a lock indicator signal when the second PLL is locked into the 
20 frequency of the input signal, and the switch may switch between the first output signal and 
the second output signal based on a status of the lock indicator signal. 

One or both of the PLLs may include a phase sensitive detector operable to receive 
the input signal and to produce a detector output signal, a loop filter operable to receive the 
detector output signal and to produce a filtered signal, and a voltage controlled oscillator 
25 operable to receive the filtered signal and to produce an oscillator signal. The phase sensitive 
detector may be further operable to receive the oscillator signal as a feedback signal of the 
PLL. It should be noted that the switch and each of the phase sensitive detector, the loop 
filter and the voltage controlled oscillator of at least one of the first and second PLLs may be 
implemented in a software process. Also, the switch and each of the phase sensitive detector, 
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the loop filter and the voltage controlled oscillator of at least one of the first and second PLLs 
may be implemented in the software process on a single digital signal processor chip. 

The phase sensitive detector, or phase sensitive detectors, each may include a Hilbert 
transformer. Moreover, the input signal may be heterodyned with a known frequency before 
5 the input signal reaches the Hilbert transformer. 

The process variable transmitter further may include an amplitude detector operable 
to sense an amplitude of the input signal and to generate a low flow signal when the 
amplitude of the input signal is below a user-controlled value. Additionally, the process 
variable transmitter may include a pre-filter operable to filter the input signal prior to 
10 processing by at least one of the first PLL and the second PLL. Based on a status of the low 
flow signal from the amplitude detector, a fixed center frequency of, for example, the second 
PLL may be switchable between the first output signal and 2nf ph , where f ph is a high cut-off 
frequency of the pre-filter. The pre-filter may be switchable between an ON state and an 
OFF state. 

15 The process variable transmitter may include a self- validating module operable to 

generate validated uncertainty parameters including a measurement value and an uncertainty 
value relating to the quality of the measurement value. The validated uncertainty parameters 
generated by the self- validating module may include a measurement status variable. The 
self-validating module may be implemented in a software process. 

20 In another implementation, a vortex flowmeter may include a flow sensor operable to 

sense pressure variations due to vortex-shedding of a fluid in a passage and to convert the 
pressure variations to a flow sensor signal in the form of an electrical signal having 
sinusoidal characteristics. This particular implementation also may include a signal 
processor operable to receive the flow sensor signal and to generate an output signal 

25 corresponding to the pressure variations due to vortex-shedding of the fluid in the passage. 
The signal processor may include PLLs having different characteristics from each other and 
operable to receive the flow sensor signal, to lock onto the flow sensor signal, and to produce 
PLL output signals indicative of the flow sensor signal. The signal processor also may 
include a switch for switching the output signal generated by the signal processor from 
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among the PLL output signals. As an example, the signal processor may be implemented by 
a software process in a digital signal processor chip. 

One of the PLLs (a "first" PLL) may be operable to lock onto the flow sensor signal 
faster than any other PLL, and another one of the PLLs (a "second" PLL) may be operable to 
lock onto the flow sensor signal with greater accuracy and greater immunity to noise than the 
first PLL. The switch may be operable to switch the output signal generated by the signal 
processor from an output signal of the first PLL to an output signal of the second PLL when 
the second PLL locks onto the flow sensor signal. When so implemented, the first PLL will 
attain lock quickly and ensure some stability during start-up and large transients, whereas the 
more accurate second PLL will ensure higher accuracy and better tracking of the vortex 
signal (and hence a better evaluation of the flowrate of the material through the passage). 

The vortex flowmeter also may include an amplitude detector operable to detect an 
amplitude of the flow sensor signal and to generate a low flow signal when the amplitude of 
the flow sensor signal is below a user-controlled value. Additionally, a filter may be 
included to filter the flow sensor signal prior to processing by the second (slower, yet more 
accurate) PLL. This filter may be switchable between an ON state and an OFF state, and is 
switched to the ON state based on the low flow signal from the amplitude detector. 

One particular method of determining a flow rate sensed by a vortex flowmeter may 
include inputting to a signal processor an input signal having sinusoidal characteristics. The 
signal processor may include a first PLL having a first bandwidth and a second PLL having a 
second bandwidth narrower than the first bandwidth. The method may include locking into 
the frequency of the input signal using the first PLL, which has a fast loop filter having a 
large natural frequency to enable the first PLL to lock quickly into the frequency of the input 
signal. Next, the method may include locking into the frequency of the input signal 
accurately using the second PLL, which has a slow loop filter having a small natural 
frequency to enable the second PLL to lock into the frequency of the input signal more 
accurately and with greater immunity to noise than the first PLL. A lock indicator signal 
may be generated when the second PLL is locked into the frequency of the input signal. 

Based the lock indicator signal, an output'of the signal processor may be switched 

from between an output signal produced by theTirst PLL and an output signal produced by 

the second PLL. Switching the output of the/Signal processor includes switching the output 
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of the signal processor from the output signal of th/first PLL to the output signal of the 
second PLL when the lock indicator signal indicates that the second PLL is locked into the 
frequency of the input signal. Additionally, thor output of the signal processor may be 
switched from the output signal of the second^LL to the output signal of the first PLL when 
the lock indicator signal indicates that the s^ond PLL is out of lock with the frequency of 
the input signal. 



The output signal of the first PLL may be provided to the second PLL as an initial 
condition frequency of the second PLL to assist lock-in by the second PLL. 

One implementation of a signal processing apparatus for acquiring a frequency of an 
10 input signal may include a first PLL having a first bandwidth and operable to lock into the 
frequency of the input signal. The signal processing apparatus also may include a second 
O PLL having a second bandwidth narrower than the first bandwidth and operable to lock into 

• : n the frequency of the input signal with greater accuracy and greater immunity to noise than 

; * the first PLL. A switch may be provided to receive signals from the first PLL and the second 

TU 15 PLL and to switch therebetween to generate an output signal. 

isy Additionally, the signal processing apparatus may include a self-validating module 

□ operable to generate validated uncertainty parameters including a measurement value 

^ corresponding to the output signal and an uncertainty value relating to the quality of the 

W measurement value. The validated uncertainty parameters generated by the self-validating 

U 20 module may include a measurement status variable. 

The multiple PLLs can be interpreted as providing a narrow band-pass filter around 
the vortex shedding frequency. The noise-rejection performance of the multiple PLLs is 
designed to give the same or better noise-rejection performance than conventional vortex 
signal processors at the minimum flow rate of the vortex flowmeter. The bandwidth of 

25 conventional vortex signal processors varies with the flow rate, giving best noise-rejection 
performance at minimum flow rate. As the flow rate increases, however, the bandwidth of 
conventional vortex signal processors increases and the noise-rejection performance is 
degraded. Note that the vortex signal has high SNR at high flow rates, so that having a 
narrow band-pass filter in high flow rates is not as important as in low flow rates, where the 

30 SNR is low. The multiple (e.g., dual) PLL approach, however, provides the same noise- 
rejection performance of conventional vortex signal processors in the minimum flow-rate (or 
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even better) over the full range of operation of the flowmeter (i.e., the bandwidth of the dual 
PLL filtering is independent of the flow rate). Therefore, some of the benefits of using dual 
PLLs instead of the conventional vortex signal processors include: better precision in high 
and low flow rates, significant improvement in the tracking performance, and extending the 
5 turn-down ratio of the vortex flowmeter. 

The details of one or more implementations are set forth in the accompanying 
drawings and the description below. Other features and advantages will be apparent from the 
description and drawings, and from the claims. 

DESCRIPTION OF DRAWINGS 

10 FIG. 1 is a diagram illustrating an exemplary structure of a process control system 

including a vortex flowmeter. 

FIG. 2 is a diagram illustrating a basic structure of the vortex flowmeter of Fig 1. 

FIGS. 3A-3C are graphs illustrating typical vortex frequency data signals. 

FIG. 4 is a diagram illustrating a basic structure of a phase-locked loop (PLL). 

15 FIGS. 5 A and 5B are a set of graphs illustrating a time-series and power spectral 

estimate of a vortex signal at a flow rate of 0. 1 1/s. 

FIG. 6 is a diagram illustrating a basic structure of a Hilbert transformer. 

FIG. 7 is a diagram illustrating one implementation of a Hilbert transformer phase 
sensitive detector. 

20 FIG. 8 is a graph illustrating spectral density of a signal after heterodyning. 

FIG. 9 is a diagram illustrating a simplified implementation of a heterodyne module. 

FIG. 10 is a graph illustrating design criteria for a Hilbert transformer receiving a 
heterodyned input signal. 

FIG. 1 1 is a graph illustrating high-pass filtering of a high frequency component of a 
25 heterodyned vortex signal corresponding to a low flow rate. 

FIG. 12 is a diagram illustrating an implementation of a vortex signal processor using 
two PLL modules. 
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FIG. 13 is a graph illustrating a start-up criterion of the vortex flowmeter design of 

Fig. 12. 

FIG. 14 is a graph illustrating start-up time vs. natural frequency co n \ of PLL1 of the 
vortex signal processor of Fig. 12. 

5 FIG. 1 5 is a diagram illustrating an implementation of a vortex signal processor using 

two PLL modules along with a prefilter. 

FIG. 16 is a diagram illustrating a system for performing center frequency analysis. 

FIG. 17 is a diagram illustrating an implementation of the switching mechanism 
associated with the vortex signal processor of Figs. 12 and 15. 

10 FIG. 18 is a diagram an implementation of the lock indicators LI1 and LI2 of Fig. 12. 

FIGS. 19A and 19B are graphs illustrating the output of the lock indicator of a PLL. 

FIG. 20 is a block diagram illustrating an implementation of a lock indicator based on 
internal signals of Hilbert Transformer Phase-Sensitive Detector. 

FIGS. 21 A and 2 IB are graphs illustrating the output of the lock indicator of Fig. 20. 

15 FIGS. 22 A and 22B are graphs illustrating the transient behavior of the lock indicator 

LI2 of the dual PLL shown in Fig. 15. 

FIG. 23 is a block diagram illustrating another implementation of a lock indicator LI2. 

FIGS. 24 A and 24B are graphs illustrating the transient behavior of the lock indicator 
LI2 shown in Fig. 23. 

20 FIG. 25 is a block diagram for comparing a single PLL followed by a smoothing 

filter, and a dual PLL structure. 

FIG. 26 is a graph comparing the performance of the single PLL followed by a 
smoothing filter, and a dual PLL structure, shown in Fig. 25. 

FIG. 27 is a flow chart illustrating a process for implementing the dual PLL shown in 

25 Fig. 15. 

FIGS. 28A and 28B are graphs illustrating outputs from the dual PLL processor 
shown in Fig. 15. 
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FIG. 29 is a flow chart illustrating a process for measuring flow of a material with a 
vortex flowmeter having the dual PLL structure shown in Fig. 15. 

FIG. 30 is a diagram illustrating features of a self- validating flowmeter 
implementation. 

Like reference symbols in the various drawings indicate like elements. 



DETAILED DESCRIPTION 

Improvements to vortex flowmeter signal processing in both digital and analog 
implementations are provided. 

Fig. 1 illustrates a simplified process control system 1 10, or flow rig, which is one 
10 implementation in which the vortex flowmeter can be used. The flow rig 1 10 may include a 
reservoir 112, a pump (e.g., a centrifugal pump) 114, conduit portions 116 (e.g., pipes having 
a 2 inch diameter), a valve 118 to control the flow rate, and a vortex flowmeter 120. Various 
other process control devices or field devices may be included in the process control system 
110, depending on the specific process control application. Furthermore, the process control 
15 system 110 may include a communication network 122 for monitoring and controlling the 
various devices of the process control system 110. 

The communication network 122 can be implemented using a typical two-wire 
communication configuration. The communication network 122 can provide analog (e.g. 4- 
20mA signals), digital, and/or pulse output communications. A central controller 124 may be 
20 provided for remote monitoring and control of the devices connected to the process control 
system 110. The central controller 124 may include a computer workstation with a graphical 
user interface for facilitating interaction between humans and the process control system 110. 

In the simplified implementation depicted in Fig. 1, the control system 110 may 
provide typical feedback control of a process. For example, the pump 1 14 may be controlled 
25 (e.g., from a signal generated at the central controller 124) to pump the material through the 
conduit portions 116 at a rate that depends on the measured flow rate sensed by the vortex 
flowmeter 120. The vortex flowmeter 120 may output to the central controller 124 a signal 
representing the flow rate of the material through the conduit portions 116. Accordingly, the 
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central controller 124 may control the pump 114 based on the programmed parameters for 
controlling and monitoring the specific process control system 110. 

Turning now to Fig. 2, an implementation of the vortex flowmeter 120 may include a 
sensor portion 210 and an electronic processor/transmitter module 212. The sensor portion 
210 may include a flow tube (e.g., a conduit-mating body) 214 for installation in coaxial 
alignment with the conduit portions 116. The sensor portion 2 1 0 may further include a bluff 
body 216 positioned to interact with the material flowing through the conduit portions 116. 
Additionally, the sensor portion 210 may include a vortex pressure sensor 218 for sensing the 
fluctuating pressure of the vortices shed by the bluff body 216 as the material flows past the 
bluff body in the conduit portions 116. This vortex pressure sensor 218 may be implemented 
using a piezoelectric sensor element or some other pressure sensitive element. 

The electronic processor/transmitter module 212 may include electronics for filtering 
and processing the signal received from the sensor portion 210. The electronic 
processor/transmitter module 212 may be implemented using one or more of, for example, a 
digital signal processor (DSP), a field-programmable gate array, an ASIC, other 
programmable logic or gate arrays, or programmable logic with a processor core. The 
electronic processor/transmitter module 212 generates a measurement of flow rate through 
the conduit portions 116 based at least on signals received from the vortex pressure sensor 
218. The electronic processor/transmitter module 212 also may include electronics for 
transmitting through the communication network 122 a data signal representing, among other 
things, the flow rate of the material flowing through the conduit portions 116. The electronic 
processor/transmitter module 212 can include electronics capable of filtering, processing, and 
transmitting digital and/or analog signals. Likewise, the central controller 124 can receive 
and process digital and/or analog signals sent from the various process control devices (e.g., 
the vortex flowmeter 120). 

In one implementation, the signal from the vortex pressure sensor 218 can be 
converted into a digital signal by an analog-to-digital converter (ADC). The digital signal 
then can be provided to a DSP that implements the electronic processor/transmitter module 
212. Accordingly, all of the processing and transmission of the vortex flowmeter can be 
accomplished digitally. Of course, the vortex flowmeter also can be implemented in an all- 
analog or mixed analog-and-digital design. 
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It should be recognized that the above-described control system 1 10 can include, for 
example, additional sensors, transmitters, power supplies, and central controllers. Moreover, 
the communication network 122 can conduct digital, analog and pulse communications 
among the various devices connected to it, and can supply power to the various connected 
5 devices (e.g., the vortex flowmeter). 

For a typical vortex sensor signal sent from the sensor portion 210, the vortex 
shedding dataXO is modeled as a sinewave of amplitude A 0 , angular frequency <j>, and with 
additive band-limited noise n(t): 

y(f) = Ao sinty, + n{t), 

10 where a>i is the frequency to be determined. With this model, it is assumed that coi oc/ the 
flow rate (though at low-flow rates the usual fluid-mechanical corrections need to be made); 

°c/ , - a well-known fluid-mechanical result - though saturation of the electronics which 
transduce and amplify the electrical signal corresponding to the pressure fluctuations might 
affect this at high flow rates; and n(t) is statistically consistent for a given flow regime, 

1 5 though it would probably increase significantly from low to high flow rates, and it is 

expected that the signal :noise ratio (SNR) improves as flow increases. The first assumption 
is common with all vortex flowmeters, the second is invoked mainly for low flow rates, and 
the third is needed to determine the resulting accuracy of the method but does not affect the 
basic design. 

20 Figs. 3A-3C illustrate some of the difficulties of measuring the vortex frequency. In 

particular, Figs. 3A-3C show vortex shedding (pressure) data recorded from a typical flow 
rig with a pipe of nominal 2 in. diameter and in which the flow is driven by impeller pumps. 

Fig. 3 A shows 0.5 seconds of data corresponding to a (medium) flow rate of 2.5 1/sec, 
acquired at a sampling rate of 1 kHz. Some conventional vortex flowmeters have been 
25 known to use peak or zero-crossing counting methods to calculate the vortex shedding 

frequency. In Fig. 3 A, the vortex frequency is about 24 Hz, which can be determined by 
counting peaks or zero-crossings. The signal shown in Fig. 3 A can be characterized as a sine 
wave with varying amplitude and frequency (note, for example, the variations in the time 
between peaks and zero-crossings). 
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The frequency estimates from zero-crossing methods have a corresponding amount of 
"jitter," that result in uncertainties in the flow measurement. However, averaging the 
frequency estimates over a number of periods may reduce this jitter to an acceptable amount. 
For example, if the vortex frequency is about 24 Hz, and flow data are to be provided every 
5 two seconds, then this will involve averaging over 48 periods of the waveform, which gives a 
reduction in variance of 48 and in standard deviation of about 7. In general, this means that 
conventional vortex flowmeters using simple zero-crossing techniques will yield acceptable 
flow data only at high flow rates. 

Fig. 3B illustrates a raw signal of a smaller flow rate of only 0.25 1/s, which is below 
10 the cut-off flow rate for most of the standard vortex flowmeters, where a standard vortex 
flowmeter is unable to measure accurately a flowrate below the cut-off point of a standard 
vortex flowmeter. It should be noted that the amplitude of the vortex shedding is much 
smaller (by a factor of 100) than the signal shown in Fig. 3 A, and that the signal is much 
noisier (e.g., from harmonic components caused by pressure pulsations from the impeller 
15 blades of the pump). A zero-crossing method would produce erroneous results due to 

spurious crossings caused by the noise. A zero-crossing method with hysteresis (such as in a 
Schmitt trigger) would remove the errors, though great care would be needed in the design of 
the hysteresis bandwidth because of the slow variations in the mean level which can cause 
the signal to move totally below or above zero (note the data around 9.5 sec in Fig. 3B). 

20 One approach is to use a low-pass (or band-pass) filter, which can lead to the data in 

the graph of Fig. 3C. The band-pass filter can be an adaptive filter using the measured vortex 
shedding frequency to establish the filter frequency settings. The basic intent of the adaptive 
filter is to establish and maintain a band-pass whose low frequency cut-off is half of the 
vortex shedding frequency and whose high cut-off is twice the vortex shedding frequency. 

25 However, failure of the zero-crossing process correctly to measure the frequency causes the 
band-pass filter to be presented with an incorrect center frequency. This induces further 
errors in frequency estimation, and lock is lost. Hence, zero-crossing processes have 
difficulty acquiring and retaining a measurement at low flow-rates. Therefore meters using 
zero-crossing methods have a minimum guaranteed flow rate. 

30 Notwithstanding that the data looks 'clean 5 in Fig. 3C, the long-term variations may 

still cause problems with a zero-crossing method. Moreover, as the frequency is now only 
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2.4 Hz (due to the lower flow rate), the method of averaging data to give accurate estimates 
will take 10 times longer (e.g., 20 seconds for the same averaging effect discussed above 
with respect to the data in Fig. 3a). The fundamental problem with a zero-crossing method is 
that it uses information at only 2 (or 3) points per period, whereas inspection of Fig. 3C 
shows that data is available over the entire period. 



PHASE-LOCKED LOOPS 

Another method of tracking sinewaves in noise uses a phase-locked loop (PLL). The 
objective of a PLL is to synchronize the frequency of a locally generated signal with that of 
an incoming signal. Data carried by the incoming signal, either in terms of frequency or 
amplitude variations, then can be accessed even if there is significant noise. The use of a 
PLL can be interpreted as narrow-band filtering around the locked frequency. Locking into 
the frequency is achieved by using a feedback loop, which uses the phase differences 
between the signals as the error. 

It is important to note that the PLL does not look at zero-crossings. Rather, the PLL 
looks for a sine wave that may be masked in the raw signal with the raw signal's spurious 
crossings caused by the noise. Essentially, the PLL transforms a raw data signal, for 
example, as shown in Fig. 3b, to a smooth data signal, for example, as shown in Fig. 3c. 

Fig. 4 illustrates an example of the basic structure of a PLL 140. The PLL 410 
includes a phase sensitive detector (PSD) 412 (also referred to as a phase-shift detector), a 
loop filter 414, a nominal center frequency input node 416, and a voltage-controlled 
oscillator (VCO) 418. The PSD 412 provides a signal e=fr#o, which is the phase error 
between the input sine wave Aosm$ t {t) and the locally generated sine wave sin$fc(f). The loop 
filter 414 is the "controller" part of the PLL feedback loop that tries to ensure that e is small 
(hopefully zero). When e is small or zero, the PLL 410 is said to be "in lock" and the 
generated sine wave s 0 (t) has the same frequency as the incoming sine wave $,-(*). 

The loop filter 414 may include a proportional-plus-integral (PI) controller 420, 

which gives a second-order closed-loop with natural frequency a) n . This natural frequency 

co n essentially determines the bandwidth of the PLL, and is used to describe the natural 

frequency of the PLL. Furthermore, the value of co n is a key design parameter: a "large" 
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value of co n yields fast tracking (e.g., during changing flow rates), but low accuracy and bad 
results for low flow rates. 

The VCO 418 generates a sine wave with the frequency proportional to the input 
voltage. The term "VCO" is merely suggestive of the analog electronic implementation. 
5 Indeed, the VCO 418 may be implemented in a software process for a digital signal processor 
(DSP). In the digital case, the VCO output is simply generated by a "sin" function in the 
DSP chip that implements the overall process. 

Finally, the signal co added in at the nominal center frequency input 416 exists to 
provide a nominal center frequency so that a sine wave is generated on start-up of the PLL 
10 and so that lock-in will be relatively quick. 

Vortex shedding frequency data have the following attributes, which are known to 
cause problems with digital and analog "standard" PLLs: 

• The vortex shedding frequency varies linearly with the flow rate. These 
variations can be as much as 100:1, depending on the size of the flowmeter, as shown 
15 in Table 1 . This is a relatively large range compared to standard PLL ranges 

(typically 1:1.5) and hence it is difficult to tune the PLL to have a small co„. 
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Nominal pipe-line size (inch)jFlow rate range (1/s) 


Frequency range (Hz) 


0.75] 0.07--2.1 


14-400 


lj 0.1--3.6 


8-300 


1.5 0.15-8.8 


3-180 


2 0.19-14.7 


2-135 


3j 0.47-32.2 


1.5-90 


4 0.89-56.7 


1-70 


6 2.4-130 


0.8-46 


8 4.6-228 


0.7-34 


10! 8.1-372 


0.6-26 


12! 12.5-540 


0.55-22 



TABLE 1: Flow rate and frequency range of the existing vortex 

flowmeter 



5 • The frequency of vortex shedding has typical values of 0.5-400Hz. This 

means a small value of a>„ is required to minimize the 2nd harmonic introduced by the 
PSD of an analog multiplier PLL, particularly at the lowest frequency of 0.5Hz, 
where G) n should be much smaller than ^ rad/s of vortex shedding. Standard PLLs 
operate most successfully at much higher frequencies (e.g., 100MHz). 

1 0 • The amplitude of the vortex shedding varies with the flow rate as A 0 =a g 2 . 

These variations are in the order of 20000: 1 . The analog multiplier PLL suffers from 
tuning problems of the PI controller, as the loop gain of the analog multiplier PLL 
depends on the amplitude of the incoming signal. 

• The vortex shedding signal is contaminated with other frequency harmonics 
1 5 generated by the pump impeller and pipe-line vibrations. The amplitude of these 

harmonics at low flow rates can be significant compared to the vortex shedding 
amplitude, as shown by the lower trace in Fig. 5, and can cause the PLL to lock into 
the wrong signal. Indeed the SNR for vortex shedding can be as low as 1 : 10 (e.g., for 
a flow of 0.1 1/s). 

20 • A PLL can lose lock and cause spurious measurements when input data are 

excessively noisy or rapidly changing. 
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These typical problems associated with vortex flowmeters have inhibited the use of analog 
and/or digital PLLs in commercial instruments such as the vortex flowmeter. Moreover, 
many standard PLLs are noise-sensitive and can cause "false-lock" (i.e., the PLLs can lock 
into harmonics of the incoming signal). 



HILBERT TRANSFORMER PHASE-SENSITIVE DETECTOR 

To avoid false-lock, give fast lock-in and ensure reliable behavior independent of 
flow amplitude, the PSD 412 may be modified to include a Hilbert transformer. In more 
detail, one of the roles of the PSD 412 in the PLL structure is to produce an error signal e(i) 
as the difference between the phase of the incoming signal and the phase of the PLL output 
signal Assuming that the incoming signal is s,{t) = A 0 s'm^ h where A 0 is its amplitude, and 
the output of the PLL 410 is s 0 (t) = sin^ 0 , then using trigonometry, the error signal - </> 0 ) 
can be produced as follows: 

x - Aosin$fios$ 0 - ^ocos$sin^ 0 = Aosm(fy - 



y = i4ocos$cos jfc + ^osin^sin^, = ^ O cos($ - 

and 



e(0 = (^-^) = arctaiJ- 



N ( sin (p i cos $ a - cos 0 . sin $ 0 ^ 



= arctan 



V 



cos (f> i cos <j> (} + sin (/> i sin <f> Q 



Accordingly, a PSD may be designed to produce an error signal - that is 
independent of the amplitude A 0 of the incoming signal. However, for this type of PSD to 
work properly, a 90° phase-shifter is required to convert A 0s m<f>i to -v4 0 cos#. While designing 
a phase-shifter for a known frequency is a relatively simple task, designing a phase-shifter to 
cover the large range of frequencies required by vortex flowmeters (e.g., 100:1) is a more 
difficult task. 

A 90° phase-shifter is simply a filter that has an output with 90° phase shift over a 
range of frequencies without affecting the amplitude of the input signal. One way to design a 
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PSD having a 90° phase-shifter is to use a Hilbert transformer. A discrete Hilbert 
transformer has a non-causal pulse response: 

7 / , 2 sin 2 0?/2) 

h( n ) = i L^ n = ±1+2,...;= 0,n = 0 

n n 

The discrete Hilbert transformer has a 90° phase shift, as its pulse response is an odd 
5 function with each element h(n)z' n having a corresponding -h(n)z n . The frequency response 
of the Hilbert transformer is: 

h(n)e- jnQ> -h(n)e jna} = -j2h(n)smnco, 

with a frequency response of the Hilbert transformer H(e f6} ) converging to -j for 0 < a>< n 
and j for n< co < 0. 

||10 To produce a practical Hilbert transformer, the frequency range is defined in terms of 

jSrj the normalized frequency, which is defined as fJf N , where /(Hz) is the frequency of the 

incoming signal and^Hz) is the Nyquist frequency (half the sampling frequency). Over the 
;.g normalized frequency range, approximations to the magnitude response of the Hilbert 

^ transformer, with unity gain, can be obtained by using, for example, the remez function in 

3 1 5 Matlab® 's signal processing toolbox. The remez function finds a best-fit (equi-ripple) to the 
□ desired magnitude response for a specified length of coefficients m and specified normalized 

% frequency range . . . /? 2 ], and then computes the corresponding parameters h(n). Note that 

the remez function tends to produce the maximum error at the transition normalized 
frequencies p\ = 0 and /? 2 = 1 . When m is even, the process produces | H(0) \ = 0 and 
20 | H(\) | = 0. When m is odd, however, the process produces | H(0) \ = 0 and \H(l)\ = 1 . 
As shown in Fig. 6, causality can be imposed by passing the incoming signal through the 
approximated Hilbert transformer 610 having length m and through a delay 612 of length 
equal to m/2, so that the phase shift between their outputs x\(t) and x 2 (t) is 90°. 

One implementation of a Hilbert transformer PSD (HTPSD) is shown in Fig. 7. The 
25 HTPSD 710 can be structured such that the input signal s,(t) is applied to the Hilbert 

transformer 712 (90° phase-shifter). The input signal sfit) also is applied to one of several 
trigonometric function nodes 714. Also input to the HTPSD 710 is the signal <f> 0 generated 
by, for example, VCO 418 of Fig. 4, where <j> 0 is transformed into components sin^ 0 and 
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cos^j. All of these signals are operated on by the various trigonometric function nodes 714 
as shown, to produce the error signal output e, which corresponds to e=# r # 0 . It should be 
recognized that the above-described HTPSD 710 can be substituted for the PSD 412 of Fig. 
4. Furthermore, the HTPSD 710 can be digitally implemented with software in a DSP. 

The HTPSD 710 shown in Fig. 7 has several advantages over standard analog and 
digital PLLs. First, the HTPSD is linear over the range of -n . . . tt. Second, the time for 
lock-in is about 8 times faster than a traditional analog PLL and about 1 .33 times faster than 
a traditional digital PLL, for all a>„. Third, the response of error e is essentially independent 
of the amplitude A 0 of the incoming signal. Fourth, a ripple-free estimate of the input 
amplitude A 0 is readily computable. Fifth, there is no second or higher harmonic component 
in the output e, when the loop is in lock. This is particularly useful in avoiding false-lock and 
in presenting a "clean" signal for the frequency estimate, but is true only if the phase-shifter 
has unity gain over the frequencies involved. 

Note that the implementation of the HTPSD 710 in Fig. 7 uses an "arctan" function 
716 to achieve proper 4-quadrant operation. Note further that the addition of a Hilbert 
transformer 712 and of the trigonometric function nodes 714 add to the computational burden 
of the DSP or other components. As such, the number of coefficients m may need to be quite 
small for fast operation. This may be accomplished if the designed frequency range; is 
restricted. Indeed, proper operation of this PSD is dependent on the Hilbert transformer 712 
design. The phase shift is always correct, but the amplitude errors caused by truncation of m 
may induce second-harmonics in the PSD output e. 

HETERODYNING 

Notwithstanding the numerous advantages of HTPSDs as described above and shown 
in Fig. 7, at least two problems may be encountered by using Hilbert transformers. First, a 
typical Hilbert transformer can cope with a frequency range of only about 1:30 with a good 
degree of accuracy while many coefficients m may be required for a large frequency range 
(recall that vortex shedding frequency range may be as much as 100: 1). Second, the 
magnitude response of the HTPSD is much more difficult to design at low frequencies, such 
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that the HTPSD becomes amplitude-dependent, which may lead to second harmonics in the 
loop. 

It is possible to overcome these problems by using heterodyning, by which the input 
frequencies are shifted to a higher range at which the HTPSD operates more effectively. For 
5 example, heterodyning the incoming signal sft) = i4 0 sin# involves multiplying sft) by 
another signal s h (t) =cos o) h t, where co h (rad/s) is a known frequency. This modulates the 
waveform to give two sidebands at co h ± to h Then, as shown in Fig. 8, a high-pass filter H h (s) 
is used to filter out the low frequency component (co h - <y,), leaving the data with frequency 
co h + Q)i for the PLL to act upon. Finally, after lock-in by the PLL, the known frequency co h is 
1 0 subtracted to obtain the measured frequency. 

Fig. 9 shows an implementation of a heterodyning module 910. As shown, the 
heterodyning module 910 includes a frequency multiplier 912 and a high pass filter 914. 
Heterodyning can be accomplished by inserting the heterodyning module 910 into the 
structure of the HTPSD 710 of Fig. 7 at a point before the signal s£t) reaches the Hilbert 
1 5 transformer 712. 

Heterodyning the signal Si(t) gives: 

A 

s, (t)s h (0 = (sin((fi>, - co h )t + <D) - sin((fi>, + co h )t + <D)), 

where O is a constant phase shift. Referring to Fig. 10, a Hilbert transformer design must be 
produced for the signal sin((<y, + 6) h )t + O) instead of sit) as discussed above with reference 
20 to Fig. 7. Again, note that the known frequency co h is subtracted from the frequency estimate 
of the PLL to obtain the frequency of interest <y,-. 

An example of a condition for selecting an appropriate heterodyning frequency co h is 
that 6) h may be chosen so that co h » co iMn and co h > o) i>max , which gives a Hilbert transformer 
design for a range of normalized frequencies: 

25 ^/,max ^;>ax + ^ 
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Therefore, with heterodyning, Hilbert transformers are designed for a smaller range of 
frequency variations with a maximum ratio of approximately 2:1, as compared to the case 
without heterodyning for which the frequency range of variations could be as much as 100:1. 

The high pass filter 914 can be designed at the minimum frequency of vortex 
5 shedding a} itmin . The low and high frequency components (co h - a>i) and (co h + are very 
close together, as also shown in Fig. 11. Accordingly, the high-pass filter 914 may be 
designed with a very sharp roll-off 1 1 10 to optimally attenuate the low frequency signal, as 
also shown in Fig. 11. For example, the cut-off frequency a> ch of the high-pass filter 914 can 
be halfway between m h - 0) ittnin and a> h + co Umin (i.e., equal to the heterodyning frequency co h ). 

1 0 Furthermore, assuming a proper design is used to implement the high-pass filter 914, 

then the heterodyning of the input signal s,{t) by the heterodyning module 910 has little or no 
effect on the SNR of the input signal Sj(t). 

MULTIPLE PLL STRUCTURE AND PROCESS 

1 5 Returning now to the two major competing issues of PLL design, i.e., fast-tracking of 

the input signal and high accuracy, it should be noted that it is difficult to improve one 
without degrading the other. In other words, the natural frequency o) n of the PLL 410 (i.e., 
the loop filter 414) produces different results depending on its value, i.e., a large co n gives 
fast-tracking, low accuracy, and the inability to stay in lock under high noise conditions, and 

20 a small co n gives high accuracy and greater immunity to noise but has a very long lock-in 
time. This problem can be overcome by using a multiple PLL structure. In one 
implementation of a multiple PLL structure, this problem, as it relates to vortex frequency 
measurement, can be overcome by using, for example, a dual PLL structure, in which each 
PLL uses a HTPSD, as shown in Fig. 12. It should be recognized that the following 

25 description of the dual PLL structure adapts quite readily to multiple PLL structures having 
more than two PLLs, with only minor modifications that will be understood by one skilled in 
the art. 
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BANDWIDTHS OF THE MULTIPLE PLL 



Fig. 12 illustrates an implementation of a dual PLL 1210, that includes a first PLL 
(PLL1) 1212, a second PLL (PLL2) 1214, and a switching mechanism 1216 for switching 
between the outputs 9\ and 02 of PLL 1 1212 and PLL2 1214, respectively. PLL1 1212 has a 
5 large natural frequency co n for fast tracking from start-up, for example to lock into the input 
signal within 1 second. In other words, PLL1 1212 may be designed to take account of start- 
up and transients by having a large bandwidth co nX . PLL2 1214, by contrast, has a small 
natural frequency co n for high accuracy and the ability to stay in lock even at low SNR. It 
should be noted that the dual PLL 1210 shown in Fig. 12 may be implemented in hardware or 
10 completely by a software process in a DSP chip. 

The center frequency of PLL1 1212 may be simply 0 O (i.e., f Q \ which is normally 

^ fixed to half way in the frequency range of the vortex flowmeter. Moreover, PLL2 1214 may 

rii be designed to have a small bandwidth cq„ 2 , to provide high accuracy (precision) in constant 

jy flow rates, which improves the uncertainty of the flow measurements. The center frequency 

:7j 15 of PLL2 1214, when out-of-lock, can be input with 0, (the frequency estimate of PLL 1 

;L 1 2 1 2), so that PLL2 1 2 1 4 is brought into lock with the help of PLL 11212. Because of the 

ffl small bandwidth a> n2 , PLL2 1214 would take a very long time to track transients with large 

fj slew-rates, and more importantly, PLL2 1214 alone might take a very long time to lock onto 

; s f the signal in start-up, which is the reason PLL1 1212 is needed. 



20 The switching mechanism 1216 can be designed to/ake the measurements (frequency 
estimate 6 X ) of PLL 1 1212 during start-up of the vortex flowmeter and during transients. 
When the flow rate is constant or varying with small j&ew-rates, however, the switching 
mechanism 1216 can take the measurements (freqifency estimate 0 2 ) of PLL2 1214 and hence 
give a more accurate estimate of the flow rate. liie switching conditions can be achieved by 
25 a careful design of a lock-indicator LI2 of PLL2 1214, so that, for example, LI2 indicates a 
signal T if PLL2 1214 is locked and '0' wh&i PLL2 1214 is out-of-lock. The lock indicator 
LI2 may be designed to be slow and siWto indicate lock-in but fast in deciding that PLL2 
1214 is out-of-lock. Moreover, once U2 indicates lock, the switching mechanism 1216 can 
switch the center frequency of PLL2 1216 to a fixed value 0 O (i.e., / 0 ). 
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The dual PLL 1210 shown in Fig. 12 may be implemented by a processor employing 
two natural frequencies 0„). The bandwidth co n] of PLL 1 1212 may be designed for fast 
tracking, whereas the bandwidth co nl of PLL2 1214 may be designed for high accuracy 
(precision). The choice of the bandwidth co n2 may be based on the precision performance of 
the vortex flowmeter 120 of Fig. 1 at the minimum flow rate. In other words, the accuracy of 
PLL2 1214 may be designed so that the signal-to-noise ratio SNR, from PLL2 1214 is 
greater than or equal to the signal-to-noise ratio SNR« of a conventional vortex flowmeter at 
a minimum flow-rate (e.g., the vortex shedding frequency corresponds to 0.5Hz). For 
example, co n2 may be at most 2 rad/s to satisfy this design criterion. 

The choice of the bandwidth o) nh however, may be based on the start-up performance 
of the vortex flowmeter 120. The design choice may assume that the pump 1 14 of Fig. 1 
starts the flow as a step, and considering the worst possible scenario, may also assume that 
the flow starts as a step to the maximum flow rate co max as shown in Fig. 13. The center 
frequency 6> 0 ofPLLl 1212 may be halfway between the maximum and minimum frequency 
15 of the vortex shedding. Hence 6 0 « co m J2. One manner of selecting the design criterion, 

therefore, is to make the start-up time t, < T (where, for example, T=l), so that the frequency 
estimate is within the lock-in range Aoj l « l.35na> nl . Accordingly: 

t - 3 Afl)2 - 3 (<p m /2-1.35«p„,) 1 , r 

4x 2 al, 4x 2 <y n 3 , 

and the bandwidth of PLL 1 1212 may be approximately: 

V' 3 

20 co , > ' nax 

In view of this, if the vortex flowmeter 120 comprises, for example, a flow tube 214 having a 
2 inch diameter (for mating with conduit portions 116 having a 2 inch diameter), then the 
maximum vortex shedding frequency is 2tc130 rad/s. Fig. 14 shows the relationship of the 
start-up time /, against different values of <y„, for this example. From Fig. 14, it can be 
25 determined that the bandwidth co M of PLL 1 1212 may be chosen so that a) nl > 20 rad/s, if 
PLL1 1212 is required to start up within 1 second, as per the present example. 
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Furthermore, as discussed above with respect to Fig. 9, after heterodyning, the known 
frequency co h may be subtracted from the frequency estimate of PLL1 1212 and/or PLL2 
1214 by, for example, a frequency subtracting module so that only the frequency of interest 
Q)i is output as the frequency estimate from PLL1 1212 and/or PLL2 1214 (depending on 
whether both or only one of the PLLs includes the heterodyning module). Also, the known 
frequency co h may be subtracted from the frequency estimate of PLL1 1212 and/or PLL2 
1214 at some point prior to inputting the signal to the switching mechanism 1216. 



PREFILTER 

As discussed briefly above, during low flow rates, the amplitude of the vortex 
shedding signal is small and the noise level becomes significant. In real applications, 
however, the noise is not simply a white Gaussian noise. There are other frequency 
harmonics generated by vibrations from, for example, the pump 1 14 and the conduit portions 
1 16, as shown in Fig. 5. Therefore, at low flow rates, a PLL may give an incorrect 
measurement by locking on to another frequency instead of the vortex signal. 

To overcome this problem, a dual PLL 1510 using^n amplitude detector 1512 to 
control a prefilter 1514 may be used, as shown in Fig. L5. The dual PLL 1510 comprises the 
PLL1 1212, the PLL2 1214, and the switching mech^iism 1216 introduced in Fig. 12. The 
prefilter 1514 may be implemented as a band-pass^filter having a high cutoff frequency^ 
(Hz) and low cutoff frequency f pl (Hz). The amplitude detector 1512 may switch on the 
prefilter 1514 if the measured amplitude A A below a threshold (user-controllable) 
amplitude A , and may switch off the prdfilter 1514 if the measured amplitude A is above a 
threshold (user-controllable) amplitude A . When the prefilter 1514 is switched on, the 
prefilter filters the input signal before the input signal it is input to the PLL2 1214. When the 
prefilter 151 4js^j^ed_off, thej^Lsjgffllis input to the PLL2 1214 unfiltered. 

The prefilter cutoff frequencies f ph and f pl may be determined using the following 
steps. First, the Reynolds number Re is used to determine the threshold flow rate Q . This is 
the flow rate below which the prefilter 1514 may be switched on. The standard relationship 
between the Reynolds number and the fluid velocity Re^VD/pi can be used, where Fis the 

fluid velocity in m/s, D is the diameter of the pipe (e.g., conduit portions 1 16) in m, and ju is 
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the fluid dynamic viscosity. Note that // is equal to 10" 6 for water. Therefore, the threshold 
flow-rate in 1/s can be expressed in terms of Re, the size of the vortex flowmeter diameter S 
in inches, and // as: 

fi=^Re. 

5 Next, the threshold amplitude A may be found using the relationship A = aQ 2 , where a is 
a calibration constant (Amplitude-Flow Ratio (AFR)), which is constant for the same vortex 
flowmeter size (diameter). Next, the high cutoff frequency f ph of the prefilter 1514 may be 
set to KQ , where K is the meter factor. The low cutoff frequency^/ may be chosen to btf p , 
= 0.5f min , where f mi „ is the minimum frequency of the vortex shedding for the vortex 
10 flowmeter 120 in use (refer to Table 1). Note that the low cutoff frequency of the prefilter 
1514 may be necessary to filter out any dc offset from any amplifier used by the electronic 
processor/transmitter module 2 12 to prepare the vortex shedding signal for measurement by 
the dual PLL 1510. 

The amplitude detector 1512 may be implemented using a peak detector to measure 
15 the amplitude A of the incoming signal. For this reason, the prefilter is switched on if A < 
A . However, regarding the peak detector, two design issues may be taken into 
consideration. First, the threshold amplitude A may be made larger than the sensor noise 
level s„. In most cases, the sensor noise is approximately 5 mV, and therefore A >5 mV. 
Second, a relay may be used to detect whether A is below A . A hysteresis of width h A can 
20 be designed so that h A = 7><j- a , which is equal to the 99.7% confidence interval of the 
amplitude estimate at the threshold level A . 

Dual PLL 1510 may be implemented such that when the prefilter 1 5 14 is switched 
on, PLL2 1214 is acting on the data (input signal) without the help of PLL 1 1212. This 
means the center frequency of PLL2 1214 may be switched to a constant center frequency 
25 6 =2nf ph , where./;,, is the high cutoff frequency of the prefilter 1514. 

As an example of a design using the above-outlined criteria, a 2 inch vortex 
flowmeter and a Re =10000 (below which the signal is very noisy) yields a threshold flow 
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rate Q =0.4 1/s, and A =0.05 V. Accordingly, the high cutoff frequency of the prefilter isf ph « 
4Hz, and the low cutoff frequency is^, « 0.5Hz. 

SWITCHING MECHANISM 

Turning now to one implementation of the switching mechanism 1216 and its design, 
5 the switching mechanism 1216 is the interface that provides the measurements of the dual 

PLL 1510, such that the center frequency of PLL2 1214 is the frequency estimate from PLL1 
1212. The switching mechanism 1216 selects between <9, from PLL1 1212 or G 2 from PLL2 
1214 as the flow measurement based on the lock indicator LI2. As discussed above, if the 
measured amplitude A is below a threshold level A , then the switching mechanism 1216 
10 output is simply 0 2 from PLL2 1214 with a fixed center frequency 6 =2nf ph , where^ is the 

high cutoff frequency of the prefilter 1514. However, if the measured amplitude A is greater 
than A , then the switching mechanism 1216 outputs the frequency estimate from PLL1 1212 
during transient and start-up, and from PLL2 1214 when the flow is constant. The switching 
mechanism 1216 selects between these modes in response to the lock indicator LI2. 




The design of the switching mechanism 12 1& may consider the following two issues. 
First, the natural frequency co ni of PLL 1 1212 ma/be selected to large enough to ensure lock- 
in of PLL2 1214 despite a noisy center frequency 0, (which is the frequency estimate from 
PLL1 1212). Second, the effect of the nois^in 0 X on the frequency estimate may be taken 
into consideration. 



20 Referring to Fig. 16, the transfer function between the center frequency G x and the 

phase noise <p 0 is 

H,(s) = ^(s)= T ' s2+T ' S - llo >^+ s 



9 X T./ +K%s + K' s 2 +2a} n2 s + o) 2 n2 ' 

where the tuning parameters are £=\, Tr2/(0 n 2 and K'=<o n2 . As shown in Fig. 16 and the 
equation above, PLL2 1214 acts on the center frequency <9, as a high-pass filter with a 
25 damping factor £ equal to 1, a natural frequency <y nl equal to a n2 , and a gain \H$a>)\ equal 
to 2/o>„ 2 for large frequencies. Assuming that co n2 is small so that \Ht(jai)\ « 2lco n2 for all 
frequencies, the variance of the phase noise due to the center frequency can be written as: 
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Var(4,) = -i-Var(ft), 

©»2 



where: 



Var(6>,) = 




Accordingly, the variance of the phase noise in PLL2 1214 due to the frequency 



5 estimate of PLL1 1212 is 



where <y„, is the bandwidth of PLL1 1212, co n2 is the bandwidth of PLL2 1214, A 0 is the 
amplitude of the incoming signal and S 0 is the spectral height of the noise corrupting the 
signal. This result demonstrates how the different parameters of the incoming signal A 0 , S 0 , 
10 and the bandwidth of PLL1 1212 and PLL2 1214, i.e., <y„, and <y„ 2 , affect the variance of the 
phase noise Var(^„) of PLL2 1214, due to a noisy center frequency 0, input to PLL2 1214. 
Assuming that the phase noise <f>„ should not exceed the capture range (7t/2) of the HTPSD 
710 used to implement the PLLs, and using the 99.7% confidence interval 3SD($,)<7i/2, the 
maximum limit of a> ni for a given a>„ 2 , S 0 and A 0 can be found as: 



As can be seen from the above discussion, the maximum limit of co nX increases as a 
result of increasing o)„ 2 . This is due to the transfer function H t (s) acting as a high-pass filter 
that rejects more noise as a n2 increases. Further, note that A refers to the threshold 
amplitude A , because below A , PLL1 1212 is disabled (i.e., switching mechanism 1216 
does not use the frequency estimate of PLL1 1212 as the flow measurement output when A 
is below A , except in extreme transient situations that cause the measured amplitude A to 
spuriously cross the threshold amplitude A ). 

Turning now to the effect of the noise in 0, on the frequency estimate the transfer 
function relating 62 to 6 X can be written as: 
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with the same tuning parameters of ^=1, Tr2/a) n2 and K'=6) n2 . As shown by the transfer 
function, //^s) is a high-pass filter with the same damping factor £ and natural frequency 
co nX of the filter Hjs) (i.e., 1 and o) n2y respectively). However the gain of the filter \HJjai)\ is 
5 unity for large frequencies and, assuming small co n2 is used, the variance of the frequency 

estimate 6*> is approximately equal to the variance of the center frequency 0, with a variance 
given by the equation provided above for Var(ft). The noise in 6 2 resulting from a noisy 
center frequency 9 X can be eliminated by fixing 0 X to a constant value, once lock is acquired 
by the PLL2 1214. As a result, the only source of noise on 9 2 is the noise in the incoming 
10 signal. 



Fig. 17 illustrates an implementation of the switching mechanism 1216 introduced in 
Figs. 12 and 15. Switches SI, S2 and S3 combine to output die frequency estimate of the 
dual PLL 1510 as well as the center frequency <9for PLL2/1214. For example, SI provides 
the frequency estimate 0 X of PLL 1 1212 as the center frequency #of PLL2 1214 when LI2 is 
15 "OFF". When LI2 is "ON," SI outputs a fixed center frequency 9 (i.e., 0 is fixed at its last 
value). S2 outputs the frequency estimate <9 2 from PLL2 1214 if LI2 is "ON" (constant flow- 
rate). Otherwise, S2 outputs the frequency estimate 0, from PLL1 1212 (e.g., during start-up 
and large slew-rate transients). S3 takes the/butput from S2 if the amplitude of the vortex 
shedding exceeds the threshold value A .Af the measured amplitude A is below A , 
20 however, S3 outputs the frequency estimate 6 2 from PLL2 1214, because PLL2 1214 has a 



^ and tl^prefil 



fixed center frequency 0=2nf ph and th£ prefilter 1512 is switched on. 



LOCK INDICATORS 

Fig. 18 shows an implementation of the lock indicators LI1 and LI2 introduced in Fig. 

12. The lock indicators LI2 may be used to ensure reliable operation of the switching 

mechanism 1216. Further, LI2 may be used to switch on a smoothing filter l/(l+sT) to filter 

the output of the dual PLL 1510 once LI2 indicates lock. Furthermore, by combining LI1 

from PLL1 1212 and LI2 from PLL2, uncertainty measurements can be generated for a self- 
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validating vortex flowmeter implementation, as discussed further below. As shown in Fig. 
18, the lock indicator 1800 may include a multiplier 1810, a low-pass filter 1812, and a relay 
1814. 

The incoming signal (PLL in put) may be s{t)=A 0 sm<frn(t), where n(t) is assumed to be a 
5 white Gaussian noise of zero mean and a spectral height S 0 . This signal passes through a 
heterodyning stage giving shAJlsinfa+nXt), where «*(/) has a spectral height of SJ2. The 
noise output from the PLL can be neglected, and the output (PLL outpu ,) is simply s„(0=sin#,. 
Accordingly, the output of the multiplier 1810 is: 

j,*(fK(0=/(0= y (cos(^,>cos(^,0)+"XO, 

10 where the spectral height of the noise n\t) is SJA. 

The signal l(t) is passed through a simple first-order filter, so that the output of the 
filter is: 

l f (0 * — cos2$ + n f (0, 

where |G(y'<y)| is the gain of the filter and is expressed as: 



15 



Accordingly, when the PLL is in lock, then $ « and the output of the filter has an 
average of A 0 /4. Heterodyning in this case has the advantage of reducing the second 
harmonic in /(f), as the frequency of the input signal is taken further by the heterodyning 
frequency co h in the magnitude response \G(ja>)\. 

20 When the PLL is out of lock, however, the filter output is: 

^(O^yl^HCcos^. -0 o )-cos(4 +&)) + «/(')» 

with a zero mean. Accordingly, the relay threshold Z* may be chosen such that the relay 
output is 1 if //>£*, and otherwise is 0. 
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In this lock indicator design, the threshold L* of the relay may be chosen to be 
A,min/4, where A 0Mn is the amplitude of the vortex shedding at the minimum flow-rate. Fig. 
19 illustrates a simulation of the lock indicator of a PLL with co=\ . The amplitude of the 
sine wave signal is A=l, and the noise added to the signal has a spectral height of le-4 and a 
bandwidth of 100 Hz. The frequency of the incoming signal is varied as a step change 
between 30 rad/s to 50 rad/s. Fig. 19 shows that Z* and the filter time constant T L should be 
chosen carefully. 



It is sufficient to use a low noise analysis^ as preftlter 1512 is used to filter out the 
noise in low flow rates, and hence high SNR is/input to PLL2 1214. Accordingly, when the 

10 PLL is out-of-lock, the noise n0) is very smzp compared to A 0 \G(jAa))\/4, and the second 

harmonic signal is reduced significantly by^the low-pass filter 1812 due to the heterodyning, 

so both noise elements can be neglected v^ien the PLL is out-of-lock. 



When out-of-lock, the amplitude of the filtered signal always should be smaller 
than the threshold of the relay 1814, so thatA 0 \G(jai)\/4 < Z*, which can be simplified to 



f A ^ 



15 1 + ^7? > ^2- 

'' I 41* 



As shown in Fig. 19, the second harmonic has a maximum amplitude around the lock- 
in frequency, where Aco L « 2o) n . As such, the filter time constant T L can be expressed as: 




If the threshold level is small and AJ4L* » 1, this can be simplified to: 



20 T L> 1 



2o)„ 



x 0 

41*, 



which shows that as L* -» 0, T L -> °o and a large time constant is needed. Therefore, this 
lock indicator depends on the amplitude of the incoming signal, and it may use a filter 1812 
with large time constant to keep the noise level, when out-of-lock, below a small threshold 
value L*. This is likely the case in the vortex filter application, as the amplitude of the signal 
25 is very small in low flow rates. 



28 



Attorney Docket No.: ^^p- 100001 
Client Ref. No.: 00,012^^ 

Another implementation of the lock indicators LI1 and LI2 is shown in Fig. 20. In 
Fig. 20, lock indicator 2000 includes trigonometric functions 2010, filters F(s), and relay 
2012. Instead of inputting the input and output of the PLL to the lock indicator 1800 as 
described above, the internal signals of the HTPSD 710 can be used to indicate lock/out-of- 
5 lock of the PLLs 1212 and 1214. Referring to Fig. 20, the signals jc=4 o sin(0 r #>) and 

y=A 0 cos(#rfa) are used as follows: the amplitude estimate A can be found using ^jc 2 + y 2 
and filtered using a simple low-pass filter F(s)=\/(l+sT L ) to obtain A F . Likewise, the signal 
y is filtered using the filter F(s) to give y F . The ratio l{f)=y r J A F can be used to indicate lock 
l{f) 1, or out-of-lock /(*)-» 0. The lock-indicator design for high SNR is acceptable for 
1 0 this application because the prefilter 1 5 1 4 is used to improve the signahnoise input to the 
dual PLL 1510 in low flow-rates. 

Fig. 21 shows the ratio l/J)—yiJ A F of a PLL with a>=l in a response to a step change 
in frequency of the incoming signal varying between 30 rad/s to 50 rad/s. The sine wave 
signal has an amplitude A^X and an additive noise with a„=0.1. 

15 Fig. 21 shows that reliable lock detection can be achieved by using a relay 2012 with 

a threshold value around 0.5 (half way between 0 and 1), so that the relay outputs "1" (i.e., 
locked) if exceeds 0.5, and the relay outputs "0" (i.e., out-of-lock) if///) is below 0.5. 
Moreover, the worst case (most noisy) of the signal 10) occurs when the PLL is out-of-lock 
and the PLL frequency estimate reaches the lock-in range, which is approximately equal to 

20 Aco L « 2co n . Therefore, the filter time constant T can be designed by considering the 

following: the signal j/=cosA^ is filtered, and the output of this filter should be kept smaller 
than the threshold level (0.5) of the relay 2012. As shown in Fig. 21, l/t) has its maximum at 
around the lock-in range Aco L . Hence: 

Jl + Aa> 2 TZ>2-+T L >^**~±, 



25 where £=1. 



An advantage of the lock indicator 2000, assuming high SNR, is that the threshold 
level L* is always 0.5 regardless of the amplitude of the incoming signal. Moreover, the lock 
indicator 2000 uses the HTPSD 710 with no second harmonic (i.e., less noise). 
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As described above, the lock indicator LI2 of PLL2 1214 may indicate out-of-lock 
during transients and start-up, so that the output of the dual PLL 1510 is simply the 
frequency estimate of PLL1 1212, and the dual PLL 1510 may output the frequency estimate 
of PLL2 1214 once the flow rate is constant and PLL2 1214 has acquired lock. 

5 The center frequency of PLL2 1214 may be the frequency estimate of PLL 1 1212, 

which is likely always locked in high flow-rates. So if the flow-rate is changing (e.g., as a 
ramp change), the lock indicator LI2 may be switching on/off during transients, as shown in 
Fig. 22. In Fig. 22, a sine wave signal is input to the PLL with its frequency varying between 
20 Hz and 50 Hz as a ramp change (rate of change of 3 1/s 2 ). The STD of the noise added to 
10 the sine wave signal is a=0.1. As shown in Fig. 22, during transients the frequency estimate 
of the dual PLL 1510 switches between the output of PLL 1 1212 and the output of PLL2 
1214, because the lock indicator LI2 of PLL2 is also switching on and off during transients. 
The reason for this is that the center frequency of PLL2 1214 is provided by PLL1 1212, 
which helps bring PLL2 1214 into lock, but PLL2 1214 soon goes out-of-lock as the flow 
jf 15 rate is still changing (for large slew-rates). 



Accordingly, a further implementation of the lock indicator (LI2) may follow the 
design of being slow and sure to indicate lock but fast in deciding that the PLL2 1214 is out 
of lock. Referring to Fig. 23, this design may be implemented by adding components to the 
lock indicator 2300. As shown, the output of LI2 (the ojdtput of the relay 1814) passes 
» 20 through a hit crossing 2310. The hit crossing 23 1 0 signals that LI2 is out of lock when the 
output of the relay changes from an ON state (i.e.yl) to an OFF state (i.e., 0). Once the hit 
crossing 23 10 indicates that LI2 is out of lock, a4witch 23 12 saves in a memory 23 14 the 



time t x when the hit crossing 2310 is ON, which is provided by the clock 23 16. A switch 

/ 

2318 works to keep the output of LI2 zero"(by outputting, for example, a "0" signal 

/ 

25 controlled by a constant module 2320) for a time t di which is a design parameter controlled 
by a time period module 2322. One/ a relational block 2324 determines that the time 
difference (determined by, for example, a multiplexer 2324 and a function block 2326) 
between the clock output t c and/ is greater than the specified value t d9 the output of LI2 may 

be simply the^>utput of the relay 1814. 



30 Adding the new components shown in Fig. 23 provides a structure that keeps the 

output of the lock indicator LI2 zero for a time t d , once PLL2 1214 is out of lock. 
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Accordingly, the dual PLL 1510 can simply output the frequency estimate of PLL1 1212 as 
the flow measurement output of the dual PLL 1510 during transients, until the time 
difference t c -t } >t d , after which the output of the dual PLL 1510 is the frequency estimate of 
PLL2 1214. The performance of the dual PLL processor is shown in Fig. 24, after adding the 
modified implementation of the lock indicator shown in Fig. 23. Note that the lock indicator 
output is OFF during transient, and not repeatedly switching between on/off as in Fig. 22. 

If the slew-rate (slope) of the flow change is not large enough to cause PLL2 1214 to 
lose lock, then LI2 would remain ON, and the dual PLL 1510 takes the measurements from 
PLL2 1214. If the slew-rate of the flow change, however, is large (maximum slew-rate is 
7T6>„ 2 72), then the dual PLL 1510 outputs the frequency estimate of PLL 1 1212 for a time t d 
until it switches to the frequency estimate of PLL2 1214. For example, t d can be chosen to be 
t d <5 9 which is an arbitrary choice. This choice means if the flow rate changes and PLL2 
1214 loses lock, then the actual output of LI2 will be taken after t d seconds. If t d is smaller 
than the transient time (and PLL2 is out-of-lock), however, then dual PLL 1510 will be 
switching between the outputs of PLL 1 1212 and PLL2 1214 during the transient, but less 
often than shown in Fig. 22. 



COMPARISON OF SINGLE PLL AND MULTIPLE PLL PERFORMANCES 

Referring to Fig. 25, the following analysis compares two PLL structures. One 
structure uses a single PLL 2510 with a bandwidth a> n] followed by a first order filter 2512 
with a time constant Tto smooth out the frequency estimate from the PLL 25 10,- and the 
other is a dual PLL 2514 structure with co nX being the bandwidth of a fast PLL (PLL1) and co n2 
being the bandwidth of a slow PLL (PLL2). The purpose of the comparison is to 
demonstrate that, with proper tuning of the dual PLL processor, the dual PLL 2514 can 
achieve superior performance over the single PLL 2510 followed by the smoothing filter 
2512. 

The transfer function between the noise ri and the frequency estimate 9 U from the 
single PLL 25 10 is: 

n' s 2 +Ks + K/T/ 
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where K=2<y„, T t =2/a>„, £,=1 and rD is a Gaussian noise of zero mean and a spectral height 
equal to SJA Q 2 , where A Q is the amplitude of the incoming sinewave signal. Now, taking the 
smoothing filter 2512 into consideration, the transfer function between ri and the smoothed 
frequency estimate 9 x /v$>\ 

H e (s) = ^(s) = 



« (s + fl)(^ +2ta„,5 + < l ) 



a<y„V 



s 3 + (a + 2oj „, )s 2 + (2a<» nl + <y 2 ).y + a<y 2 , ' 

where a=\/T. 

The variance of the frequency estimate d^is given by: 



4 ' 



1 0 which may be simplified to yield the following expression for the variance of the smoothed 
frequency estimate: 

Var(0 ]f )= f'^ 2 S 0 . 

This expression can be used to explore the effect of the smoothing filter 25 12 on the 
frequency estimate of the single PLL 2510: 

15 (1) If a=oo (i.e., T=0 and no filtering is used), then 

Var( ^ ) ^r^' s o =Var( ^> <Xiy »'- 

4^4 0 

(2) If a large time constant T is used, and a « co„ u then 

(3) If a=co nU then 

\6A 0 4 
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and a reduction by a factor of 4 is obtained, compared to the "no filtering" case above. 

Assuming that the accuracy (of the software process or hardware) of the dual PLL 
2514 is equal to the accuracy (of the software process or hardware) of the single PLL 2510 
followed by the smoothing filter 2512, such that Var(6^)=Var(6'i / ) ? then: 

f 1/^2 V /3 



i/r 



For example, for a 2 inch diameter vortex flowmeter, co nX (of, e.g., PLLl or single 
PLL 25 10) may be greater than or equal to 20 rad/s to achieve a start-up within I second, as 
discussed above with respect to Fig. 14. Fig. 26 shows the relationship between co n2 and the 
required time constant T of the smoothing filter 2512 so that Var(6^)=Var(6 f 1/ ) 5 for a given 

1 0 <£>„,=20 rad/s. Fig. 26 further indicates that for single PLL 25 10 (which can also apply to 

PLLl in the dual PLL 2514) having a bandwidth a) nl =20 rad/s, a smoothing filter 2512 with a 
time constant of approximately 7M.5 seconds may be required to smooth out Q x to get the 
same variance of the frequency estimate obtained using the dual PLL 2514 with a> n2 =2 rad/s. 
Note that a larger oo nX would require a larger time constant T to smooth out the data and keep 

1 5 the variance of the frequency estimate equal to that obtained using the dual PLL 2514 with 
a) n2 =2. Having a smaller co nX (e.g., 15 rad/s), however, would require a smaller time constant 
7 1 , but this degrades the start-up performance, as shown in Fig. 14. 

Considering the transient performance of the single PLL 2510 followed by the 
smoothing filter 2512, assume that a frequency ramp with a slope a, is input to both the 
20 single PLL 2510 structure and the dual PLL 2514 structure. For the single PLL 2510 

structure, the total error is the combined error due to the single PLL 2510 and the smoothing 
filter 2512. The single PLL 2510 is shown to have an error equal to 2aJo) nX and the error 
resulting from a first order filter F{s)= U(l+sT) due to a ramp input is equal to: 

UmsE(s) = lim— sT =a f T . 
■ v ->° s \ + sT ' 

25 Therefore, the total error of the single PLL 25 1 0 structure following a ramp frequency 

is: 
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The dual PLL 2514 structure, however, has an error equal to 2a/co nl , provided a 
proper design of LI2 is achieved (i.e., a small filtering time constant is used in LI2). 
Therefore, the transient performance of the single PLL 2510 followed by the smoothing filter 
5 25 12 has an additional error during transient equal to oijT. Thus, for the same accuracy 

performance of both the single PLL 2510 (plus the smoothing filter 2512) and the dual PLL 
2514, the dual PLL 2514 gives better start-up and transient performance than the single PLL 
2510 (with smoothing filter 2512), particularly with a large slew-rate a f . 



1 0 A PROCESS FOR DESIGNING THE MULTIPLE PLL 

Fig. 27 illustrates one implementation of a process 2700 for designing, for example, 
the dual PLL 1510 (shown in Fig. 15). First, an appropriate bandwidth for the PLL1 1212 
(i.e., the "faster" PLL) may be chosen (step 2710). As discussed above, the PLL1 1212 may 
be designed to lock into the vortex frequency quickly, (e.g., within 1 second). Accordingly, 
15 the natural frequency co nX may be chosen according to: 



where co max is the maximum frequency of the vortex shedding signal. 

Next, the appropriate bandwidth for the PLL2 1214 (i.e., the "slower" PLL) may be 
chosen (step 2712). As discussed above, the PLL2 1214 may be designed to lock into the 
20 vortex frequency with better noise-rejection performance than conventional vortex 

flowmeters. This can be achieved, for example, by selecting co n2 such that co n2 < 0A5co min , 

where o) mm is the minimum frequency of the vortex shedding signal (which is, from Table 1, 
approximately 2Hz for a 2 inch size flowmeter). 

Next, the lock indicator LI2 may be designed so that, for example, the output of the 
25 dual PLL 1510 may be simply the output of PLL 1 1212 during start-up, and transients with 
large slew-rate (step 2714). Once PLL2 1214 acquires lock, the output of the dual PLL 1510 
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may be simply the output of PLL2 1214 when the flow rate is constant or during flow 
variations with small slew-rates. The lock indicator based on the internal signals of the 
HTPSD 710 (as discussed with respect to Fig. 20) may be used, and the time constant of the 
filter may be based on T>l/a>„ 2 . 

5 Next, the bandwidth of the prefilter 1514 may be designed so that it filters out the 

frequency harmonics generated by the impeller pump 1 14 and other interference at low flow 
rates (step 2716). In most cases, these harmonics have high frequencies and may be outside 
the bandwidth of the prefilter 1514. The choice of the high frequency cut-off ^„ (Hz) may be 
based on: 

10 /, A =K^%SRe, 

where K is the meter ^-factor, ju is the dynamic viscosity, S is the size of the vortex 
flowmeter (diameter) in inches and Re is the Reynolds number below which the prefilter 
1514 is switched on. The low cut-off frequency of the prefilter 1514, however, may be 
simply^, « 0.5f min where f min is the minimum frequency of the vortex shedding signal in Hz. 

1 5 Finally, the relay (e.g., amplitude detector 1512) that switches on the prefilter 1514 

may have a hysteresis having a width h r to take account of the amplitude noise (step 2718). 
The hysteresis width may be chosen such that h r = 3a- , that is the 99.7% confidence 

interval of the amplitude noise at the threshold level A , which can be estimated during 
calibration. 

20 Outputs from a dual PLL processor implemented according to the above discussions 

are shown in upper and lower traces of Fig. 28. As shown in the upper trace of Fig. 28, the 
output of the dual PLL is the flow-estimate from PLL1 during start-up, and the output of the 
dual PLL switches to PLL2, once locked, and remains locked during transients with small 
slew-rate. As shown in the lower trace of Fig. 28, in transients with larger slew-rates, the 

25 output of the dual PLL is from PLL1 until PLL2 acquires lock, and then the output of the 
dual PLL is simply taken from PLL2. 
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A PROCESS FOR MEASURING FLOW RATES USING THE VORTEX 
FLOWMETER 

Referring to Fig. 29, an exemplary flow rate measurement procedure 2900 
implemented by the vortex flowmeter 120 may be used to measure the flow rate of material 
5 through the conduit portions 116. In general, the procedure may be performed continuously 
to provide real-time data to the central controller 124 of the process control system 110 
shown in Fig. 1 . 

During start-up of the process control system 110 5 material flow may be induced in 
the conduit portions 1 16 by the pump 1 14 (step 2910). At this time, the sensor portion 210 of 

10 the vortex flowmeter 120 may begin to sense the flow of the material (i.e., sense vortices 
being shed in the wake of the bluff body 216 as the material passes) and to send electrical 
signals to the electronic processor/transmitter portion 212, where measurement of the vortex 
shedding frequency will commence (step 2912). Before or after the raw vortex shedding 
signal reaches the electronic processor/transmitter portion 212, it may be amplified to 

1 5 facilitate processing. Further, the raw vortex shedding signal may be converted to a digital 
signal by an analog-to-digital converter (ADC) before or after it reaches the electronic 
processor/transmitter portion 212. In any case, the vortex shedding signal may be input to 
the electronic processor/transmitter portion 212 for further processing. 



20 (vortex shedding signal) as an input, the amplitude detector 1512 may determine whether the 



25 flow prefilter 1514 may be controlled to be "OFF" (step 2916). 

In this normal (higher) flow situation, PLL1 1212 will be used to lock quickly into the 
vortex shedding frequency (steps 2918 and 2920). Once PLL1 1212 has locked into the 
vortex shedding frequency, its estimated output f\ may be output as the flow measurement 
signal of the vortex flowmeter 120 (at least until PLL2 1214 can lock into the vortex 
30 shedding frequency more accurately) (step 2922). 



Once the electronic processor/transmitter portion 212 receives the electrical signal 



amplitude A of the incoming signal is above a predetermined (user-controlled) threshold 
amplitude A (step 2914). As discussed above, a low amplitude may indicate a low flow 
situation, and a higher amplitude may indicate a normal (higher) flow situation. If the 



measured amplitude A is higher than the predetermined threshold level ( A ), then the low 
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It may be beneficial to set the center frequency / of PLL2 1214 equal to the 
frequency estimate/, of PLL1 1212 to "jump-start" the lock-in of PLL2 1214 once PLL1 
1212 is locked in (step 2924). Further, PLL2 1214 may be actively locking in while PLL1 
1212 is locking in, if the dual PLL 1510 is designed to have both PLL1 1212 and PLL2 1214 
5 operational and existing concurrently. Next, PLL2 1214 may be checked to see if it is in lock 
(steps 2926 and 2928). For example, lock indicator LI2 can be used to indicate whether 
PLL2 1214 is in lock, as discussed above with respect to Figure 12. If LI2 is "ON," then 
PLL2 1214 is locked into the vortex shedding frequency input to the electronic 
processor/transmitter portion 212. Accordingly, if LI2 is "ON," the center frequency / of 
10 PLL2 1214 may be controlled to remain at its last value (step 2930), and the output of the 
vortex flowmeter 120 may be set to be the frequency estimate f 2 of PLL2 1214 (step 2932). 

Generally, during start-up, the flow of the material will ramp up quickly so as to 
cause the measured amplitude A of the vortex shedding signal to be greater than the 
threshold amplitude A of the amplitude detector 1512. As such, during start-up and during 
15 large transients, the exemplary flow rate measurement procedure 2900 may apply steps 2914 
through 2932, which make sequential use of both PLL1 1212 and PLL2 1214 to arrive at the 
flow measurement output. 

In low flow rate conditions, it may be appropriate to measure the flow rate using the 
more accurate, but slower, PLL2 1214. In particular, if the measured amplitude A (Step 

20 2914) of the incoming signal is below the predetermined (user-controlled) threshold 

amplitude A , then the low flow prefilter 1514 is controlled to be "ON" (step 2934). At this 
point, the flow rate is determined to be a low flow rate, and the most appropriate procedure 
for accurately measuring the flow rate may be to use PLL2 1214 exclusively. Accordingly, 
the center frequency / of the PLL2 1214 is set to the fixed center frequency f Q of the dual 

25 PLL 1510 (step 2936). Further, the estimated output/ 2 of PLL2 1214 is output as the flow 
measurement signal of the vortex flowmeter 120 (even when PLL2 1214 has not yet 
accurately locked into the vortex shedding frequency) (step 2922). PLL2 1214 is 
continuously checked to see if it is in lock (steps 2926 and 2928). Again, if LI2 is "ON," 
then PLL2 1214 is locked into the vortex shedding frequency input to the electronic 

30 processor/transmitter portion 212. Accordingly, if LI2 is "ON," the center frequency / of 
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PLL2 1214 may be controlled to remain at its last value (step 2930), and the output of the 
vortex flowmeter 120 may be set to be the frequency estimate^ of PLL2 1214 (step 2932). 

It should be noted that, for example, steps 2914 through 2936 may be repeated, 
selectively skipped, rearranged, grouped, nested looped, set into subroutines or randomized 
5 as appropriate. Furthermore, the steps of the exemplary flow rate measurement procedure 
2900 may be supplemented with various other steps for measuring the flow rate of the 
material through the conduit portions of the process control system 110 with the vortex 
flowmeter 120. 

It should be recognized that while the implementation discussed above relates to a 

10 vortex flowmeter having two PLLs, one of which locks faster than the other, the vortex 

flowmeter could instead include any number of PLLs having different characteristics (e.g., 
locking speed) with respect to each other. Furthermore, in the case of implementing the 
PLLs in a DSP, for example, using PLL software processes, it is possible to use a single PLL 
software process and then merely adjust the software parameters used by the PLL software 

1 5 process to implement each of various faster or slower PLLs. In other words, the PLL 

software process in a DSP might use a first set of parameters (e.g., values) to implement a 
fast-lock PLL, estimate the vortex frequency using the fast-lock-PLL-tuned PLL software 
process, output from the vortex flowmeter the vortex frequency estimated by the fast-lock 
PLL, tune the PLL software process to implement a slower-lock PLL by changing the 

20 parameters used by the PLL software process, use the fast-lock PLL frequency estimate as a 
starting point center frequency of the slower-lock PLL, and estimate the vortex frequency 
using the slower-lock-PLL-tuned PLL software process. Once the slower-lock PLL is locked 
onto the vortex frequency, the DSP can output from the vortex flowmeter the vortex 
frequency estimated by the slower-lock PLL. Additionally, instead of using a switching 

25 software routine or device that switches between the outputs of two actual PLLs 

(implemented by software or hardware), the "switching" may take the form of changing the 
parameters used by the PLL software process. Moreover, this process of tuning the PLL 
software process can be repeated as necessary to implement any number of PLLs having 
different characteristics. 

30 Furthermore, many of the steps in the exemplary processes shown by Figs. 27 and 29 

can be rearranged, supplemented with other steps, combined or selectively removed. 
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SELF-VALIDATING SENSOR ATTRIBUTES 



The vortex flowmeter may use in a contcbl system that includes self- validating 
sensors. To this end, the vortex flowmeter may be implemented as a self-validating meter. 
5 Self-validating meters and other sensors are/escribed in U.S. Patent No. 5,570,300, entitled 
"SELF-VALIDATING SENSORS", whicK is incorpo^ 

In general, a self- validating meter provides, based on all information available to the 
meter, a best estimate of the value of a parameter (e.g., flow rate) being monitored. Because 
the best estimate is based, in part, on nonmeasurement data, the best estimate does not always 
10 conform to the value indicated by the current, possibly faulty, measurement data. A self- 
validating meter also provides information about the uncertainty and reliability of the best 
estimate, as well as information about the operational status of the sensor. Uncertainty 
information is derived from known uncertainty analyses and is provided even in the absence 
of faults. 

15 A self- validating meter may provide four basic parameters: a validated measurement 

value (VMV), a validated uncertainty (VU), an indication (MV status) of the status under 
which the measurement was generated, and a device status. The VMV is the meter's best 
estimate of the value of a measured parameter. The VU and the MV status are associated 
with the VMV. The meter produces a separate VMV, VU and MV status for each 
20 measurement. The device status indicates the operational status of the meter. 

The meter also may provide other information. For example, upon a request from a 
control system (e.g., a request from the central controller 124), the meter may provide 
detailed diagnostic information about the status of the meter. Also, when a measurement has 
exceeded, or is about to exceed, a predetermined limit, the meter can send an alarm signal to 
25 the control system. Different alarm levels can be used to indicate the severity with which the 
measurement has deviated from the predetermined value. 

VMV and VU are numeric values. For example, VMV could be a temperature 
measurement valued at 200 degrees and VU, the uncertainty of VMV, could be 9 degrees. In 
this case, there is a high probability (typically 95%) that the actual temperature being 
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measured falls within an envelope around VMV and designated by VU (i.e., from 191 
degrees to 209 degrees). 

In one implementation, the electronic processor/transmitter module of the meter 
generates VMV based on underlying data from the sensors. First, the electronic 
processor/transmitter module derives a raw measurement value (RMV) that is based on the 
signals from the sensor. In general, when the electronic processor/transmitter module detects 
no abnormalities, the electronic processor/transmitter module has nominal confidence in the 
RMV and sets the VMV equal to the RMV. When the electronic processor/transmitter 
module detects an abnormality in the sensor, the electronic processor/transmitter module 
does not set the VMV equal to the RMV. Instead, the electronic processor/transmitter 
module sets the VMV equal to a value that the electronic processor/transmitter module 
considers to be a better estimate than the RMV of the actual parameter. 

The electronic processor/transmitter module generates the VU based on a raw 
uncertainty signal (RU) that is the result of a dynamic uncertainty analysis of the RMV. The 
electronic processor/transmitter module performs this uncertainty analysis during each 
sampling period. Uncertainty analysis, originally described in "Describing Uncertainties in 
Single Sample Experiments," S.J. Kline & F.A. McClintock, Mech. Eng. , 75, 3-8 (1953), has 
been widely applied and has achieved the status of an international standard for calibration. 
Essentially, an uncertainty analysis provides an indication of the "quality" of a measurement. 
Every measurement has an associated error, which, of course, is unknown. However, a 
reasonable limit on that error can often be expressed by a single uncertainty number 
(ANSI/ASME PTC 19.1-1985 Part 1, Measurement Uncertainty: Instruments and 
Apparatus). 

As described by Kline & McClintock, for any observed measurement M, the 
uncertainty in M, w Mi can be defined as follows: 

where Mis true (M irue ) with a certain level of confidence (typically 95%). This uncertainty is 
readily expressed in a relative form as a proportion of the measurement (i.e. w h4 /M). 

In general, the VU has a non-zero value even under ideal conditions (i.e., a faultless 
sensor operating in a controlled, laboratory environment). This is because the measurement 
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produced by a sensor is never completely certain and there is always some potential for error. 
As with the VMV, when the electronic processor/transmitter module detects no 
abnormalities, the electronic processor/transmitter module sets the VU equal to the RU. 
When the electronic processor/transmitter module detects a fault that only partially affects 
the reliability of the RMV, the electronic processor/transmitter module typically performs a 
new uncertainty analysis that accounts for effects of the fault and sets the VU equal to the 
results of this analysis. The electronic processor/transmitter module sets the VU to a value 
based on past performance when the electronic processor/transmitter module determines that 
the RMV bears no relation to the actual measured value. 

To ensure that the control system uses the VMV and the VU properly, the MV status 
provides information about how they were calculated. The electronic processor/transmitter 
module produces the VMV and the VU under all conditions-even when the sensors are 
inoperative. The control system needs to know whether VMV and VU are based on "live" or 
historical data. For example, if the control system were using VMV and VU in feedback 
control and the sensors were inoperative, the control system would need to know that VMV 
and VU were based on past performance. 

The MV status is based on the expected persistence of any abnormal condition and on 
the confidence of the electronic processor/transmitter module in the RMV. The four primary 
states for MV status are generated according to Table 2. 



Expected 


Confidence 


MV Status 


Persistence 


in RMV 




not applicable 


nominal 


CLEAR 


not applicable 


reduced 


BLURRED 


Short 


zero 


DAZZLED 


Long 


zero 


BLIND 



Table 2: Primary states of MV status 
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A CLEAR MV status occurs when RMV is within a normal range for given process 



conditions. A DAZZLED MV status indicates that RMV is quite abnormal, but the 
abnormality is expected to be of short duration. Typically, the electronic 
processor/transmitter module sets the MV status to DAZZLED when there is a sudden 
change in the signal from one of the sensors and the electronic processor/transmitter module 
is unable to clearly establish whether this change is due to an as yet undiagnosed sensor fault 
or to an abrupt change in the variable being measured. A BLURRED MV status indicates 
that the RMV is abnormal but reasonably related to the parameter being measured. For 
example, the electronic processor/transmitter module may set the MV status to BLURRED 
when the RMV is a noisy signal. A BLIND MV status indicates that the RMV is completely 
unreliable and that the fault is expected to persist. 

Two additional states for the MV status are UNVALIDATED and SECURE. The 
MV status is UNVALIDATED when the electronic processor/transmitter module is not 
performing validation of VMV. MV status is SECURE when VMV is generated from 
redundant measurements in which the electronic processor/transmitter module has nominal 
confidence. 

The device status is a generic, discrete value summarizing the health of the meter. It 
is used primarily by fault detection and maintenance routines of the control system. 
Typically, the device status is in one of six states, each of which indicates a different 
operational status for the meter. These states are: GOOD, TESTING, SUSPECT, 
IMPAIRED, BAD, or CRITICAL. A GOOD device status means that the meter is in 
nominal condition. A TESTING device status means that the meter is performing a self 
check, and that this self check may be responsible for any temporary reduction in 
measurement quality. A SUSPECT device status means that the meter has produced an 
abnormal response, but the electronic processor/transmitter module has no detailed fault 
diagnosis. An IMPAIRED device status means that the meter is suffering from a diagnosed 
fault that has a minor impact on performance. A BAD device status means that the meter has 
seriously malfunctioned and maintenance is required. Finally, a CRITICAL device status 
means that the meter has malfunctioned to the extent that the meter may cause (or have 
caused) a hazard such as a leak, fire, or explosion. 
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One way in which to implement the vortex flowmeter as a self-validating meter is to 
exploit other signals generated by the vortex flowmeter in addition to the frequency 
measurement and its associated uncertainty, such as the lock-indicator of PLL1 (LI1), and the 
lock-indicator of PLL2 (LI2). LI1 and LI2 are combined to give self-validating 
measurements and metrics VMV and VU measurements. In one implementation, the flow 
rate measurement Q is found using Q = Kf h where K L/s/Hz is the calibrating factor, andy; 
Hz is the frequency of the vortex shedding. The flow rate measurement Q is associated with 
an uncertainty AQ so that the true measurement Q true lies in the interval 

Q~kQ<Qtrue<Q + &Q 

with a certain degree of confidence (95%, for example). This interval can be expressed in 
terms of a relative uncertainty measure 

w Q = &Q/Q as 0(1 - w Q ) < Q true < Q{\ + w Q ). 

Assuming K and f t are independent, then the uncertainty in the flow rate measurement 
of the vortex flowmeter can be found using a propagation rule, which gives the following 
result: 



AQ 
Q 



SK 



AK 
Q) 



SQ 

ft 



Yf 



4£ 
Q 



1/2 



K 



f, 



which shows that the uncertainty in the flow measurement Q is the sum of the uncertainties 
in the calibrating K-factor K and the frequency estimate f h where the uncertainty of/J can be 
found by estimating the variance of the frequency from the vortex flowmeter, on-line (95% 
20 confidence interval is equal to approximately 2oy). The flowmeter calibrating factor and its 
uncertainty (±0.5% to ±2%) are normally calibrated by the manufacturers, and hence it is 
invaluable to have a high accuracy process measuring the frequency of the vortex shedding 
as it reduces the measurement uncertainty of the flow rate. 



25 



43 



Attorney Docket No.!^B2- 100001 
Client Ref. No.: 00,0 \T 



MV status 


LI1 


LI2 


Comments 


CLEAR 


ON 


ON 


High-flow rate 


BLURRED 


OFF 


ON 


Low flow rate 


BLURRED 


ON 


OFF 


Transients 


DAZZLED 


OFF 


OFF 


Searching 


BLIND 






If DAZZLED stays longer than / seconds 



Table 3: MV metrics for the vortex flowmeter 



5 Table 3 summarizes an implementation of MV metrics for the self- validating vortex 

flowmeter shown in Fig. 30. CLEAR indicates that the VMV measurement is good, and the 
O flow measurements are in the high flow rate region, where the accuracy is expected to be 

3 high. The CLEAR signal is generated if LI 1 of PLL 1 is "ON" and LI2 of PLL2 is "ON." 

12 BLURRED indicates that the VMV measurement is good, and the flow measurements is 

f U 10 within the low flow rate region or in a transient condition. This signal is generated, for 
! t y example, when the prefilter 15 14 is switched on, and hence when LI1 is "OFF" and LI2 is 

!L "ON" (low flow rate) or when LI1 is "ON" and LI2 is "OFF" (transient or start-up). 

2 DAZZLED indicates that the VMV measurement is in a search state, and is generated if LI1 

IjlJ and LI2 are both "OFF." BLIND indicates that the VMV measurement is not good, and the 

15 flowmeter is unable to provide any measurements due to, for example, a fault in the 
electronics (e.g., the power supply is off), or caused by no flow passing through the 
flowmeter. It is simply generated when the DAZZLED state stays for longer than / seconds, 
where t is a time parameter fixed by the designer. 

In another implementation, a single lock indicator signal (e.g., LI2) could be used to 
20 generate self-validated uncertainty parameters for the vortex flowmeter. For example, if the 
lock indicator is ON, then the self-validating meter can indicate a CLEAR MV status, and if 
the lock indicator is OFF, then the self-validating meter can indicate a DAZZLED MV status. 
Further, if the lock indicator remains OFF for the predetermined wait-for-recovery amount of 
time. 

25 A number of implementations have been described. Nevertheless, it will be 

understood that various modifications may be made. For example, the dual PLL 1510 may 
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be implemented by a combination of hardware (e.g., circuitry) and software (e.g., as a 
software process in a DSP chip), on an element-by-element basis. As another example, the 
vortex flowmeter electronic processor/transmitter portion 212 may be implemented having a 
multiple PLL structure, i.e., more than two PLLs. 
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